Idealized glass transitions for a system of dumbbell molecules 
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The mode-coupling theory for ideal glass transitions in simple systems is generalized to a theory for 
the glassy dynamics of molecular liquids using the density fluctuations of the sites of the molecule's 
constituent atoms as the basic structure variables. The theory is applied to calculate the liquid-glass 
phase diagram and the form factors for the arrested structure of a system of symmetric dumbbells 
of fused hard spheres. The static structure factors, which enter the equations of motion as input, 
are calculated as function of the packing fraction tp and the molecule's elongation £ within the 
reference-interaction-site-model and Percus-Yevick theories. The critical packing fraction <p c for 
the glass transition is obtained as non-monotone function of ( with a maximum near ( = 0.43. 
A transition line is calculated separating a small-£-glass phase with ergodic dipole motion from a 
large-^-glass phase where also the reorientational motion is arrested. The Debye- Waller factors at 
the transition are found to be somewhat larger for sufficiently elongated systems than those for 
the simple hard-sphere system, but the wave-number dependence of the glass-form factors is quite 
similar. The dipole reorientations for £ > 0.6 are arrested as strongly as density fluctuations with 
wave vectors at the position of the first sharp diffraction peak. 

PACS numbers: 64.70.Pf, 61.20.Lc, 61.25.Em 



I. INTRODUCTION 



The mode-coupling theory (MCT) for idealized liquid- 
glass transitions has been proposed originally as a micro- 
scopic approximation theory for the dynamics of simple 
liquids JjJ. The MCT equations formulate the idea that 
correlation functions for density fluctuations have to be 
evaluated self consistently with the correlation functions 
for force fluctuations. The derived equations require the 
static structure factor as input, which is anticipated to 
be a smooth function of the wave vector and of con- 
trol parameters like the packing fraction cp. The equa- 
tions exhibit a bifurcation singularity for certain values 
of the control parameters, say, for if — ip c . For ip < ip c , 
the solutions describe ergodic liquid dynamics, while for 
tp > ip c nonergodic dynamics is obtained describing an 
amorphous solid. The arrested glass structure for ip > tp c 
is characterized by glass form factors, also referred to as 
nonergodicity parameters. They generalize the concept 
of the Edwards- Anderson parameter from the theory of 
spin glasses [|J; they can be determined in scattering 
experiments and molecular-dynamics-simulation studies. 
The MCT equations can be solved by asymptotic expan- 
sion using, e.g., \ip — ip c \ as a small parameter || |j. The 
leading order results establish universal results for the 
glassy dynamics ||. Anticipating that these universal 
formulas are valid also for mixtures and for molecular 
glass-forming systems, extensive tests of the MCT with 
data from experiments and simulations have been carried 
out during the past ten years ]^| ([). Due to the inven- 
tion of improved spectrometers and progress in simula- 
tion techniques, the work of testing MCT is still an active 
field. Let us mention as particularly impressive recent 
examples the studies with the optical Kerr effect] |, [|, 
the depolarized- light-scattering work for toluene EC | , the 
quantitative tests of the form factors for silica and 



the scaling-law analysis of simulation data for a polymer 
model |12| . The indicated tests suggest that MCT deals 
properly with some essential features of glass-forming 
systems. 

There is a problem in the tests of the universal MCT 
formulas: the range of validity of these leading order 
asymptotic results is not universal. For example, the 
time interval for the density-fluctuation decay according 
to von Schweidler's power law depends non-trivially on 
the wave vector of the fluctuations. Fitting data by a 
power law for times outside the regime of validity of the 
asymptotic law may be possible but can yield misleading 
conclusions |l3| . The range of validity of the leading or- 
der result can be determined by calculating the leading 
corrections or by comparing with the numerical solutions 
of the full equations of motion || [I| . But for such dis- 
cussions one has to analyze the complete equations of 
motion, i.e. one needs an understanding of the micro- 
scopic details of the system. Thus, there is the necessity 
to extend the MCT so that models can be analyzed which 
describe the experimental situation closely. This is the 
motivation for the present paper where MCT shall be 
extended to molecular liquids and where this extension 
shall be exemplified for a hard-dumbbell system. 

Extensions of MCT to molecular systems have been 
studied already, generalizing the concept of a density- 
fluctuation correlator to the one of infinite matrices of 
correlation functions formed with tensor-density fluctu- 
ations 0, H [l§ H £| g || . The re- 
sults calculated for the glass-form factors for a model 
of water |u| |l8| and for a liquid of linear molecules j^] 
could be used to explain simulation data quantitatively. 
Promising results for anomalous oscillation spectra for 
a dipolar- hard-sphere system have been calculated ]po[ . 
For the model of a dilute solute of linear molecules in a 
solvent of spherical particles, the MCT equations could 
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be fully solved El], g2|. The solutions were used to 
demonstrate the applicability of the universal formulas 
also for reorientational motion and to explain the charac- 
teristic difference between the a-peaks for dielectric-loss 
and depolarizcd-light-scattering spectra, as they have 
been observed in some experiments for van-der-Waals liq- 
uids. The cited work shows that MCT studies may con- 
tribute to the understanding of glassy dynamics which is 
beyond the implications of universality. 

The MCT equations based on the tensor-density de- 
scription of molecular systems have a different mathe- 
matical structure than the ones studied so far. It is un- 
clear whether the bifurcation dynamics of these equations 
exhibits the same universal laws as derived within the 
MCT for atomic systems. It is not obvious that codes 
can be developed for the numerical solutions of these 
equations within the regime of glassy dynamics. There- 
fore it was suggested to base the MCT for molecular sys- 
tems on the site- representation |24|, |25| . This leads to 
equations with n-by-n matrices where n is the number of 
atoms producing the force centers in the molecules. For 
the simple case of a dilute solution of linear molecules, 
it was shown that this approach yields results pij| in 
semiquantitative agreement with the much more involved 
tensor-density theory |^2). In the present paper, this 
work shall be continued with the intention to demon- 
strate a complete set of results for the glassy dynamics 
of a system of linear molecules. 

The paper is organized as follows. The basic general 
MCT equations are obtained in Appendix [a| by modifi- 
cation and generalization of the previous work ||| p4[ . 
They are specialized in Sec. II to a formulation of the 
equations of motion for the coherent and incoherent 
density correlation functions for the symmetric-hard- 
dumbbell system. The static structure factors, which 
determine the mode-coupling coefficients, are evaluated 
within the RISM theory. To analyze their features in 
Sec. HI, they are decomposed in their various angular 
momentum contributions which are evaluated within the 
Percus-Yevick theory. Section |y] explains the phase di- 
agram for the system and the glass-form factors. The 
conclusions (Sec. |v|) summarize the results while the dis- 
cussion of the correlation functions is left for a following 
paper p6|. 



II. A MODE-COUPLING THEORY FOR A 
SYSTEM OF SYMMETRIC DUMBBELLS 



Vector e, = (r- 



? )/L abbreviates the axis of the 



ith molecule. Denoting the mass of atom a as m a , the 
total mass M — ttia + ms and the moment of iner- 
tia / = mAiriBL 2 /M determine the thermal velocities 
Vt = ^/ksT/M and vr = ^/ksT/I for the molecule's 
translation and rotation, respectively. Here T denotes 
the temperature. Let us introduce also the center-of- 
mass position rf = (m^f^ 4 -I- ms^ B )/Af and the co- 
ordinates z a of the atoms along the molecule's axis: 
za = L(itib/M) 7 zb = —L(rriA/M). The basic struc- 
tural variables are the two interaction-site-density fluc- 
tuations for wave vectors q: 



N 



P% = ex P W • rt), a = A or B. 



(1) 



The site-site static structure factors S^ b 



(pTp\)I n 



provide the simplest information on the equilibrium 
structure of the system. Here (• • • ) denotes canonical 
averaging. Because of isotropy, S^ b depends only on the 
wave number q = \q\. The site-site static structure factor 
Sq h consists of the intramolecular and the intermolecular 
parts. The former is denoted as w^ b ; for a rigid dumbbell 
molecule it is given by Wg b — S ab + (1 — S ab ) j (qL). Here 
and in the following ji(x) denotes the spherical Bessel 
function of index i. The static structure factors S^ b 
shall be combined to a two- by- two matrix S g , and sim- 
ilar matrix notation will be used for other site-site cor- 



relation functions, 
tion M, & S Q 



The site-site Ornstein-Zernike equa- 

- r,,r-l _ nr. 1~ 1 Qdb 



PCq 



relates to the 



site-site direct correlation function c ab . 

The structural dynamics of the system shall be de- 
scribed by the interaction-site-density correlators 



K b (t) = ( P a Atyp»)/ N - 



(2) 



These are real even functions of time obeying F^ b (t) — 
F ba (t). The short-time expansion can be written as 



Fg(*) 



In 2 ] 
2 1 J « 



t 2 + o(t 4 



(3) 



The continuity equation reads p°~ 



iq ■ ji, where 



the longitudinal current fluctuation is given by j'| = 
J2i exp(iq ■ r°-) with v? denoting the velocity of atom 
a in the ith molecule. Therefore, one gets J° b = ((<f • 
i,f )*(<? ' j'5))/iV(7 2 , whose explicit expressions for a rigid 
dumbbell molecule are given by [E9| 



A. The model 

A system of N rigid dumbbell molecules distributed 
with density p is considered. The molecule shall be 
described within the interaction-site formalism [^7], p8[ , 
where the constituent atoms shall be called atoms A and 
B. Let r.°-, a = A or B, denote the position vectors of 
the atoms in the ith molecule, so that L = Ir^ 4 — rf\ 
denotes the distance between the two interaction sites. 



Jf = v 2 T wf+v R (lz a z b ) {5 ab +(l~6 ab ) Uo(QL)+MqL))]. 

(4) 

The dynamics of the tagged molecule shall also be con- 
sidered. It is described by the self part of the interaction- 
site-density correlators 



O) = </>r,.(*)Vr,.(o)>. 



(•5) 



Here p°j s = exp(iq ■ r s a ) with r s a denoting the position 
vector of atom a in the tagged molecule. The short-time 
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expansion of the correlator ¥ QjS (t) is given by Eq. (|) with 
S g replaced by w q . The same function J q b determines 
the short-time dynamics of F q h s (t) since the velocities 
of different molecules at the same time are statistically 
independent. 

For later convenience, it shall be shown here how the 
correlation functions in the interaction-site representa- 
tion can be expressed in terms of the ones in the tensor- 
density description. Following the convention in Refs. |^ 
and ^2], coherent tensor-density fluctuations p™(q) for 
the angular-momentum index £ and the helicity index m 
shall be defined by decomposing the ith molecule's po- 
sition variable in plane waves exp(iq ■ ff) for the center 
of mass ff and in spherical harmonics Y™(ei) for the 
orientation vector e): 

N 

P T(q) = i e >/4^exp(tf. ff ) (6) 

i=l 

The structural dynamics is described by the matrix of 
correlators 

*2,(?,t) - (pT(M*p?(.M)/N\ go = (o,o, 9 ). 

(7) 

The general correlators {p™(q,t)*p™ (q, 0)) can be writ- 
ten as linear combination of the functions $>%g,(q, t); they 
vanish for m 7^ m! if q = <fo |2l|] . In particular, the equi- 
librium structure is described by the static correlation 
functions 

SrAq) = (pT(Qo)*P i P(qo))/N. (8) 

Since the position vectors of the interaction sites can be 
written as ?® = rf + z a ei, the Rayleigh expansion of the 
exponential in Eq. (|l|) yields the formula 

p% = ]T V2t+lMqz a ) p° e (qo)- (9) 

1 

Substitution of this expression into Eq. (0) leads to an 
expression for the density correlators in the site represen- 
tation in terms of those in the tensor-density description: 

Ff(t) = J2 v / (2^+l)(2^' + l) U (qz a ) n ,(qz b ) $° u ,{q,t). 

1,1' 

(10) 

In particular, the site-site static structure factors S" are 
related to the tensorial ones via 

Sf = VW+W* + T}Mqz a )jt>(qZb) S%,{q). 

(11) 

Similarly, one obtains formulas relating tagged-molecule 
correlators in the site representation and those in the 
tensor-density description: 

F£{t) = ^v/(2£+l)(2£' + l) 

xje(qz a )ji>(qz b )$° s>u ,(q,t), (12) 



where <&® u ,(q,t) denotes the self part of <I>° £ , (q,t). Since 
^° s ,u'(l^) = $w, one gets 

wf = + l)ji(qz a ) j t (qz b ). (13) 

e 

In the present paper, a system of symmetric dumb- 
bell molecules, consisting of two fused hard spheres of 
diameters cLa = ds = d and masses tua = tub = M/2, 
shall be considered. The elongation parameter £ = L/d 
quantifies the bond length. All equilibrium properties of 
such a hard-dumbbell system with a fixed elongation are 
specified by the packing fraction 

<P = pVo, 1 ^ = 5^ 3 ( 1 + |C-|C 3 ). (14) 

Here Vq is the volume of a dumbbell molecule. Through- 
out the rest of this paper, the diameter of the constituent 
atom is chosen as the unit of length, d = 1. 

For a symmetric system, there are only two indepen- 
dent density correlators, since F q A {t) — F q B {t). It is 
convenient to perform an orthogonal transformation to 
fluctuations of total number density p^ and "charge" 
density p%: 

p%=(pf±p*)/V2, x = NorZ. (15a) 
The transformation matrix P = P _1 reads 

p = 7lQ-I)' <15b) 

It diagonalizes the matrices S g , w q , and J q : 

(P S, P) xv = S xy S x = Sf A ±S AB , (15c) 
(Pw,Pr=P<, tflf = l±j («0, (15d) 
(P3 q P) xy = 5 xy {v 2 T w x q 

+ \v 2 R e [lT(7o(«?C)+Ja(gC))]}, (15e) 

where x, y = N or Z. Also the matrix of density correla- 
tors is diagonalized. Introducing the density correlators 
4> q (t) normalized to <f>®(t = 0) = 1, one gets 

r q (t)={pWpm)/ Ns ^ ( i6a ) 

(PF g (t)P)-¥ = **»^(t)S* (16b) 

and similar equations hold for the normalized tagged- 
molecule correlators 4> q s (t): 

^, s W = «>WVf, s (0))/<, (17a) 

(PF 5 , s (t)pr=f^,(f)^. (17b) 

There is an additional property due to the symmetry 
of the molecule. Since the intermolecular parts of S AA 
and S AB are the same, one gets S? — w q . A similar 
reasoning for the charge-density correlators leads to 

<(*)=<C(*)- (18) 

For a system of symmetric dumbbell molecules, there is 
only one independent coherent density correlator, viz. 
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B. MCT equations for the density correlators 

The MCT equations of motion for the density correla- 
tors consist of an exact Zwanzig-Mori equation and the 
approximate expression for the relaxation kernel in terms 
of the mode-coupling functional, whose derivation is de- 
scribed in Appendix |A|. For a system of symmetric dumb- 
bells, these equati ons can be simplified considerably pi] . 
Multiplying E qs. (|A3|) -(|A7| ) from left and right with P 
given by Eq. (15b) and inserting 1 = PP between every 
pair of matrices, all equations are transformed to diago- 
nal ones. Thus, there are two sets of equations, one for 
<f)q (t) and another for <t>^ (i) . As explained in connec- 
tion with Eq. (|l^), the charge-density correlator 4>g(t) is 
identical to its self part, (j>^ s (t), which shall be treated 
separately below. Thus, the only correlator describing 
the coherent density fluctuations is the total density cor- 
relator (f>^ (i), whose Zwanzig-Mori equation reads 



d t 2 <(*) + (^) 2 <(t) 

+ (nf) 2 / dt'm»(t-t')d t ,4>%(t') = 0.(l9 a ) 
Jo 

The characteristic frequency 51^, which specifies the ini- 
tial decay of the correlator by <j>% (t) = 1 - |(0^ tf + 
0(t )j is given by 

«) 2 = <7 2 {4 [i+io(0] 

+ y R C 2 ll-Jo(q()-.n(q()}}/S?. (19b) 
The relaxation kernel reads (t) = [<j) N (t)\, where 



Eqs. (A£)-(A7) lead to 



1 



F?\J\ = 2j dk VN{ ^ kJ) ~ h ^ (2 ° a) 
V N (n- h n 1 - __P_QN oN oN 



x{q-{kc% (20b) 



with p = q — k and — 2c AA . One gets from Eq. ( All ) 
for the nonergodicity parameters = <fi^ (t — > oo): 



f?=F?[f N }/{l+F?[f N ]} 



(21) 



Notice that Eqs. (19a), (|20|), and ( |2l| ) are formally identi- 
cal to the corresponding equations for simple systems || : 
the difference is in the definition of the correlators and 
of the direct correlation functions. In particular, one can 
show that the preceding equations ( 19 )— (|2l|) reduce to 
the ones for simple systems in both the £ — > and £ — > oo 
limits. 

The MCT model for the hard-dumbbell system 
(HDS) will be defined by two further technical as- 
sumptions. First, the site-site structure factors S® 
and the direct correlation functions c ab are evalu- 
ated within the reference-interaction-site-model (RISM) 
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FIG. 1: Phase diagram of the symmetric-hard-dumbbell sys- 
tem where the packing fraction is denoted by ip and the elon- 
gation parameter by (. The solid curve marks the type-B 
liquid-glass-transition line, <p c = ¥>c(C)- The dashed curve 
denotes the type-A-transition line between phases II and III, 
<Pa = <pa (C)- The type-A-transition line terminates at the 
critical elongation £ c = 0.345 marked by an arrow. The hor- 
izontal arrow marks the transition point of the hard-sphere 
system (HSS). 



integral-equation theory (27|, £8| [30[. Second, the wave 
numbers are discretized to 100 equally spaced values 
q = 0.2, 0.6, 1.0, • • • , 39.8. The details of the transfor- 
mation of the mode-coupling functional to a polynomial 
in the discretized variables can be found in Ref. [| 

The discussion of Eq. (|2l]) can follow that considered 
previously for simple systems ||. For a given £, one finds 
a critical packing fraction tp c = ip c (() so that = for 
ip < ip c and > for tp > ip c . Figure |l| exhibits 
the control-parameter plane for our system; the full line 
represents the cp c -versus-C curve. The regime I, i.e., the 
states with (£, tp) below the full line, are the liquid states. 
For states on and above the line, the density-fluctuation 
dynamics is nonergodic. It is the purpose of this paper 
to explain the origin of this liquid-glass-transition curve 
and to quantify the arrested glass structure. 

Comments on some technical details of our calculations 
might be in order. We solved a set of equations in the 
RISM theory to obtain using the non-equally spaced 
wave-number grids introduced in Ref. |H} The resulting 
Sq b has been subsequently transformed to the one on the 
above-mentioned equispaced grids using a cubic spline in- 
terpolation J32]| . Our results for the HDS are based on the 
Sg b so obtained. Occasionally, we will refer to results for 
the hard-sphere system (HSS), i.e. the HDS with £ = 0. 
For consistency, calculations for the HSS have also been 
done using a static structure factor which is based on 
the same numerical method just mentioned. However, 
the numerically obtained structure factor for the HSS is, 
due to the interpolation procedure, slightly different from 
the one of the analytic Percus-Yevick theory p8| . Since 
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the transition is sensitively dependent on the structure- 
factor peak, this leads to slightly different results for the 
HSS from the previous ones reported in Ref. ||, where 
the analytic Percus-Yevick theory is used. Typically, the 
differences are less than 1%. Therefore, the results for 
the HSS based on the two different static inputs can 
be regarded as essentially the same. The nonstandard 
wave-number grids from Ref. |3l] and the subsequent in- 
terpolation procedure have been adopted because of the 
following reason. The present model shall be extended to 
one where the constituent atoms carry opposite electrical 
charges. Thereby it will be possible to study the inter- 
play of steric-hindrance effects and Coulomb-interaction 
effects. The method developed in Ref. [IT] is well suited 
for treating such a system in which both the short- and 
long-ranged interactions are simultaneously present. To 
have the results for the HDS as reference model, it seems 
adequate to carry out the calculation of the static input 
function strictly within the same frame. 



MCT equations for the tagged-molecule 
correlators 



One gets the Zwanzig-Mori equation for the normal- 
ized tagged- molecule correlator (j) x ^ s {t) (x = N or Z) 
for a symmetric-hard-dumbbell system by transforming 
Eq. (AS) as explained above for deriving Eq. (19a): 



(t) becomes identi- 



diX iS (i-i')9t'C s (*')=().(22a) 



The characteristic frequency CI* s specifies the initial de- 
cay of the correlator by c/) x s (t) = 1 - \{Cl x jf + 0(i 4 ), 
and it is given by 



W = i 2 K+XC [lTjo(9C)Tj 2 (gC)] / [i±jo(gC)]}. 

(22b) 

The relaxation kernel can be written as m x s (t) — 
F x [(j) x (t),(j) N (t)}, where Eqs. @ and (|ATo|) lead to 



r« , = _F q 

9 ' sL/s ' iJ 16^ 3 q 2 



dkl q ' P \ fJf^cN. 



( c p ) w kSp fk,sfp> 

(23) 



with p = q — k. From the long-time limits of Eqs. ( |A8| ) 
and (A9), one gets for the nonergodicity parameters 



f«, s =^Af?>f}/{ 1 +Ks[f?>r}}- ( 24 ) 

The mathematical structure of the two sets of the equa- 
tions (^2|) (|24|) for x = N and Z is the same as the one 
studied previously for the tagged-particle-density corre- 
lator in a simple liquid In particular, the set of equa- 
tions for x = N reduce to the one for the tagged-particle- 
density correlator in both the £ — > and £ — > oo limits. 



In the same limits, the correlator . 
cally zero. 

Equations for the tagged symmetric dumb- 

bell immersed in a liquid of symmetric dumbbells are also 
formally identical to the ones treated in Refs. |24| and |25| 
for the symmetric dumbbell molecule dissolved in a sim- 
ple liquid. This is because the coherent density fluctua- 
tions of the surroundings in the former case is character- 
ized only by the correlator 4> q (t), i.e., a scalar correlator, 
and this feature is shared with the latter case. In analogy 
to the findings in the previous studies, one finds a line of 
transition points (fA — ¥?a(C)) provided £ < £ c = 0.345. 
This line, which is shown in dashed in Fig. separates 
glass states in regime II and regime III. In regime II, the 
reorientational motion is ergodic, i.e., the states deal with 
the amorphous analogue of a plastic crystal. In regime 
III, also the reorientational motion is nonergodic, since 
fqs > 0. Crossing the dashed line by, e.g., increasing </?, 
fq S change continuously (type- A transition). Crossing 
the heavy line, f q changes discontinuously (type-B tran- 
sition). The interest of the present studies concerns the 
transition from the liquid to a glass with all density corre- 
lators arrested, as obtained for £ > Q c by increasing tp. As 
a representative situation with strong steric hindrance for 
reorientational motion, molecules with £ = 1.0 shall be 
analyzed in detail. For £ approaching Q from above, the 
steric hindrance for reorientations weakens and molecules 
with £ = 0.4 shall be used to demonstrate this case. 



III. STRUCTURE OF THE RELAXED SYSTEM 



The static structure factor for the total density fluctu- 
ations, , is the basic input of our theory. It quantifies 
the simplest information on the averaged particle distri- 
bution, anticipating the system to be relaxed in a canon- 
ical equilibrium state. The latter is assumed to be an 
amorphous one. It may be metastable, e.g., with respect 
to crystallization. The variation of with changes of 
the packing fraction <p and the molecule's elongation £ 
provides the key for explaining the phase diagram in 
Fig. |l|. Extending earlier work (33| to the high density 
regime, S q shall be analyzed in this section. 



A. Static structure factors and angular correlations 

Figure || exhibits results for S q calculated from the 
RISM theory (27[ |2^, ^] for the two representative elon- 
gations ( = 0.4 and 1.0 at and near the critical pack- 
ing fraction <£> C (C)- ^ ov small q, the structure factor is 
small. Because of the dense packing, the compressibility 
K q cx Sq for long-wavelength fluctuations is strongly sup- 
pressed. These fluctuations are irrelevant for the glassy 
arrest in our system. The phase diagram does not change 
more than 1% if fluctuations with, say, q < 3 are cut off. 
Fluctuations with, say, q > 10 are relevant, since — 1 is 
of order unity. But in this regime, the structure factor is 
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FIG. 2: Static structure factor Sq for the total density 
fluctuations as function of wave number q for the elongations 
£ == 0.4 (upper panel) and 1.0 (lower panel). The results refer 
to packing fractions ip — <p c (l + e) with e — (solid lines), 
e = -10~ 5/3 (dashed lines), and e = +10" 5/3 (dotted lines). 
Here ip c denotes the critical packing fraction; it is given by 
tp c = 0.675 and 0.565 for £ = 0.4 and 1.0, respectively. The 
first sharp diffraction peak in Sq for £ = 0.4 occurs at q = 7.0 
in the discretized wave-number grids, and its heights are 4.54, 
5.47, and 6.61 with increasing (p. The corresponding peak for 
£ = 1.0 occurs at q = 7.4, and its heights are 3.33, 3.75, and 
4.24 with increasing ip. Here and in the following figures the 
diameter of the spheres is used as unit of length, d = 1. 



not very sensitive with respect to changes of the density. 
Therefore, the liquid-glass transition is driven mainly by 
the changes of Sg for q w 7, i.e. by the fluctuations with 
wave vectors near the position of the first sharp diffrac- 
tion peak. This feature is analogous to the one found in 
the hard-sphere system [Q. However, the results shown 
in Fig. ^| are for molecular systems in which angular cor- 
relations should play an important role as well. This 
subsection is devoted to discuss how angular correlations 
manifest themselves in . 

To proceed, let us decompose the in terms of the 
spherical-harmonic expansion coefficients S^, (q) defined 
in Eq. (||): the coefficient 5*Q (g) describes the static 
center-of-mass density fluctuations, and the higher co- 
efficients probe the angular correlations. This decom- 
position can be derived from Eq. (Jll|) by noticing the 



definition S, 



N 



s, 



AA 



s, 



AB. 



S? = 2V(2£+1)(2£' + 1) 



f f :cvcn 



xji(q(/2) Je ,(qC/2)S u ,(q). (25) 
Here the angular-momentum indices £ and £' take only 
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FIG. 3: Spherical-harmonic expansion coefficients Su(q) of 
the structure factor for the elongations £ = 0.4 (upper panel) 
and 1.0 (lower panel) at the critical packing fraction p = <p c 
for the angular-momentum indices I = 0, 2, and 4. 



even numbers due to the top-down symmetry of the 
dumbbell molecule. It is clear that the coefficients S® t , (q) 
contain more information than since the latter can 
be expressed in terms of the former, but not vice versa. 

The expansion coefficients Sf$, (q) have been calculated 
within the Percus-Yevick (PY) theory |54|, |3f| up to the 
angular-momentum-index cutoff £ C ut = 6. For the sym- 
metric dumbbell, this results in 30 independent coeffi- 
cients to be dealt with in solving the PY equation. The 
representative results at the critical packing fraction for 
the diagonal coefficients S^(q) are shown in Fig. |[ For 
the small elongation £ = 0.4, the density fluctuations for 
q w 7 are dominated by those of the center-of-mass de- 
grees of freedom, £ — 0, while contributions from the 
reorientational correlations are rather small. On the 
other hand, for the large elongation £ = 1.0, the static 
structure for q « 7 is primarily caused by the reorien- 
tational function S^iq), while the center-of-mass com- 
ponent Sqo(?) on ly shows a weak structure. A strong 
peak at q w is also seen in the coefficient S^iq) for 
C = 1.0, which is a precursor of a nematic instability. 
The increased importance of the higher coefficients for 
larger elongations is demonstrated even more clearly by 
comparing S^(q) for the two elongations. These features 
of the coefficients S® e (q) for small and large elongations 
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FIG. 4: The solid lines denote the static structure factor 
S^ 1 for the total density fluctuations for the elongation £ = 
0.4 (upper panel) and £ = 1.0 (lower panel) at the critical 
packing fraction ip = ip c (0- The dashed and dotted lines 
are the results of decomposition based on Eq. ( pB[ ) , where the 
numbers in parentheses (1,1) indicate the component in the 
decomposition (see text). 



are in accord with those found in Ref . |3(| albeit for fluids 
of hard ellipsoids in which the aspect ratio plays a role 
similar to (1 + £). 

Figure ^ exhibits the decomposition of S^ at the crit- 
ical packing fraction based on Eq. (Eq). The solid lines 



denote Sq calculated from the RISM theory, and the 
dashed and dotted lines denote the terms in the decom- 
position using the coefficients S^(q) from the PY theory. 
Cross terms (£ ^ £') are omitted to avoid overcrowding of 
the figures. The function from the RISM theory and 
the one based on Eq. ( p5|) with the coefficients S" e ,(q) 
from the PY theory are found to be in good agreement 
with each other, and therefore it makes sense to discuss 
the decomposition of using the results from two dif- 
ferent integral-equation theories. As mentioned in con- 
nection with Fig. |2|, the glass transition of our system is 
driven by the first peak in S^ centered at q 7. Figure ^ 
shows that, for the small elongation, the first peak is pri- 
marily determined by the center-of-mass density fluctua- 
tions, whilst the contributions from higher order angular 
correlations are responsible only for the peaks located in 
the higher-g region. On the other hand, when the elon- 
gation is large, the contribution from the center-of-mass 
degrees of freedom gets suppressed, but the higher or- 
der angular correlations become much more important 
for determining the first peak: the first peak is primar- 
ily accounted for by the (2, 2) contribution. Thus, the 



static density fluctuations determining the cage for the 
glass transition are of different origins for small and large 
elongations, respectively. 

A comment shall be added concerning the strong peak 
in S^iq) at q « shown in Fig. H for £ = 1.0. It is clear 
from Eq. ( f25|) that the function S% 2 (q) contributes to S^ 
with a prefactor j2(9C/2) 2 , which is proportional to q* 
for small q. As a result, this strong peak in S^iq) hardly 
contributes to Sj^ in the small wave-number regime; it 
only gives rise to a small peak centered at q ss 3.5 as 
shown in Fig. [|. Also, it is seen that the (0, 0) component 
has a peak at the same q range. However, it is found 
that the small peak at q 3.5 is canceled out by the 
(0, 2) component, which is not shown in the figure. All 
this together results in the small and flat Sff for q < 5. 
We conclude that the strong peak in S^iq) at q w 
is irrelevant for the glass formation for the elongation 
£ = 1.0 within our theory. 



The intramolecular correlation functions 



(x 



N,Z) from Eq. (15d) enter the mode-coupling vertices 
implicitly via the site-site Ornstein-Zernike equation for 

"231 



S„ and explicitly via Eq. 



be decomposed in analogy to Eq 



Using Eq. (13), they can 
5|): 



,N (Z) _ 



I : even (odd) 



i)M«C/2) 5 



(26) 



Wq starts for q — at the value 2 (0) and then it os- 
cillates for q > 5 around the value 1. The first oscillation 
minima of occur at q = 11 A and 27.4 for ( = 0.4, 
and at q — 4.5 and 11.0 for £ = 1.0. For Wq , the first 
minima are located at q = 19.4 and 35.0 for £ = 0.4, and 
at q = 7.8 and 14.2 for C = 1.0. 

Let us consider the change of the structure factor as 
function of the elongation £ for fixed packing fraction 
ip. Figure || exhibits the result for if = 0.56. It is 
seen that for small elongations (the upper panel) the 
first peak height decreases with increasing the elonga- 
tion, whilst the opposite trend is seen for large elonga- 
tions (the lower panel). This feature can be explained 
in terms of the spherical-harmonic expansion coefficients 
S( e (q) as follows. As discussed in connection with Fig. |J, 
the center-of-mass density fluctuations (I = 0) are pri- 
marily responsible for determining the first peak in 
for small elongations, whereas the angular correlations of 
the index i = 2 are more relevant for large elongations. 
The strength of the center-of-mass correlation becomes 
weaker as the elongation is increased, and this explains 
the decrease of the first peak height in for small 
elongations. On the other hand, when the elongation 
is large, £ = 2 component is relevant, and this angular 
correlation becomes stronger with increasing the elonga- 
tion. This explains the increase of the first peak height 
in with increasing £ for large elongations. Thus, the 
non-monotonic ^-dependence of the first peak height in 
Sq for the fixed packing fraction shown in Fig. || is due 
to the different origin of that peak for small and large 
elongations. 
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FIG. 5: Static structure factors for the total density 
fluctuations at the fixed packing fraction ip — 0.56 for various 
elongations as indicated in the figure. 



B. Preferred orientations for nearest neighbors 

A digression might be adequate for a better under- 
standing of the equilibrium structure of our molecular 
systems, especially, of angular correlations for nearest 
neighbors. Such angular correlations can best be in- 
vestigated through the molecular pair correlation func- 
tion g(Vi2, 0i, 02, ^12)- Here ri2 denotes the center-to- 
center separation, and the three angles 0i, 02, and ^12 

1 — (f> 2 specify the relative orientations of the two linear 
molecules in the so-called r-frame. The pair correlation 
function can be expanded as |35| 



3(^12,01,02,012) 



4tt 



E 



xr^(0i,0i)r;;(0 2 ,0 2 )*.(27) 

The <?™£ 2 (ri2) are the spherical harmonic expansion co- 
efficients, and can be calculated within the PY the- 
ory 35 1 . In the present work, the coefficients are 
calculated up to cutoff £ cut — 6. 

The center-to-center radial distribution functions 
3oo( r i2) fo r representative elongations at the large pack- 
ing fraction ip — 0.56 are shown in Fig. |6[ along with the 
radial distribution function for hard spheres (C = 0.0) at 
the same packing fraction. As the elongation increases, 
the first peak position in <7oo( r i2) increases, the height 
of the peak decreases, and the peak becomes broader 
and somewhat irregular, with a shoulder developing at 
separations just beyond 7*12 = 1. For the elongation 
C = 1.0, the shoulder turns into a broad prepeak cen- 
tered at r\2 ~ 1.1. 



FIG. 6: Center-of-mass component of the molecular pair 
correlation function, <?oo(?~i2), as function of the center-to- 
center distance Tyz at the fixed packing fraction ip — 0.56 for 
various elongations 



It becomes more difficult to interpret the g™i 2 {r 12) for 
non-zero values of l\ and £2- Therefore it seems ade- 
quate to follow Streett and Tildesley |3j| and consider 
cuts through the space of the four variables determin- 
ing the function g in Eq. (27). Typical cuts for dis- 
cussing the relative orientations of two linear molecules 



are 
tt/2," 
(0i = 
(0i = 
orientation 



12 — 
% = 

02 = 



[33): (i) the "T-shaped" orientation (fit = 0, 6 2 = 
any value), (ii) the "crossed" orientation 
/>12 = t/2), (iii) the "parallel" orientation 
tt/2, 4>i2 = 0), and (iv) the "end-to-end" 
9i = 2 = 0, 0i2 = any value). These ori- 
entations lead to efficient packing at close approach in 
the sense that they all lead to the close contact of the 
constituent atoms, and thus contribute to the first peak 
in . Most of the orientations at high densities can be 
broadly classified as being similar to one of these four. 
Because of the computational reason to be described in 
the next paragraph, the "crossed" and "parallel" orien- 
tations shall be combined to define the "CP-type" orien- 
tation (0i = 02 = 7r/2) by averaging over the angle 4>i 2 : 



3(^12,01 



1 

2^ 



(28) 



Notice that the corresponding 0i2-averaged pair corre- 
lation functions for the "T-shaped" (0i = 0, 2 = tt/2) 
and "end-to-end" (0i = 02 = 0) orientations remain the 
same as the original ones since the angle <p\2 is irrele- 
vant in defining these two orientations. Figure |?] exhibits 
representative results. 
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FIG. 7: The (^12-averaged molecular pair correlation function 
g(ri2, 01,82) as denned in Eq. ( Bgh for the "T-shaped" (Si = 0, 
#2 = 7r/2) and the "CP-type" (9i = 62 = 7r/2) orientations 
(see text) at the packing fraction tp — 0.56 for the elongations 
C = 0.4 (upper panel) and 1.0 (lower panel). 



Before embarking on the conclusions to be drawn from 
Figs. H and (?], let us make comments on the cutoff prob- 
lem in the summation in Eq. (|27|). To check the con- 
vergence, we also performed the same calculations with 
4ut = 4. It is found that for distances m > 1 + (, 
the series converges rapidly, i.e., the difference between 
the results with £ cut = 4 and those with £ cut = 6 is small. 
However, for distances 7*12 < 1 + the difference is rather 
large reflecting the slow convergence of the series, and 
this effect becomes more pronounced with increasing den- 
sity and elongation. An indication of the lack of conver- 
gence is that the functions for £ = 1.0 shown in the lower 
panel of Fig. ^ can take unphysical negative values in the 
range rn < 1 + C The reason for the slow convergence in 
that range is that the function g[rx2,0t,6%, $12) has step- 
like features because of the hard-core repulsion, which 
cannot be accurately represented by a truncated series 
with a small number for £ CU f This problem is likely to 
be a feature of the spherical harmonic expansion for any 
model in which the molecule has a relatively hard asym- 
metric core. It is also found that a better convergence 
is achieved for the 0i2-averaged correlation function de- 
fined in Eq. ( p8| ) than the original g(r 12, 9%, #2, 4>i2), an< 4 
this is why we have chosen the averaged ones to display 
the results. Despite these unwelcome features, it is antici- 
pated that qualitative features of the angular correlations 
are captured even with £ cut = 6. Notice that the men- 
tioned cutoff problem do not influence the results to be 
presented for the MCT since those are based solely on 
Sjf calculated from the RISM theory. 

The increase of the most probable nearest neighbor 



center-to-center separation with increasing elongation, 
which is demonstrated in Fig. ^|, suggests that the major- 
ity of nearest neighbor pairs adopt orientations for which 
the center-to-center distance of closest approach increases 
with increasing elongation. Therefore, it seems likely 
that "CP-type" orientations do not contribute heavily 
to this peak in <7oo(?~i2), because their closest approach 
remains in the region ri2 ~ 1 irrespective of the elonga- 
tion. It is also clear that "end-to-end" orientations are 
unimportant, because their minimum approach distance 
(ri2 = 1 + C) lies well beyond the distance at which the 
first maximum occurs in ffoo( r i2)- Hence, the major con- 
tributions to the first peak in <7oo( r i2) are likely to come 
from orientations of the "T-shaped" type and ones close 
to it. The first maximum in the g(r 12, 8 'i, #2) for the "T- 
shaped" orientation for £ = 0.2, 0.4, 0.6, 0.8, and 1.0 oc- 
cur at r\2 — 1.13, 1.24, 1.33, 1.41, and 1.46, respectively, 
as shown in Fig . |t] for £ = 0.4 and 1.0. These positions are 
very close to the first maximum positions in <7oo( r i2) for 
each elongation shown in Fig. ^|. This evidence is consis- 
tent with a strong predominance of "T-shaped" nearest 
neighbor orientations. We therefore conclude that at the 
most probable nearest neighbor distance there is a strong 
preference for "T-shaped" orientations over all others. 

We next consider how the "CP-type" correlations man- 
ifest themselves in the 500(^12)- As can be inferred from 
Fig. fj], such correlations would lead to a peak centered at 
r\2 « 1 irrespective of the elongation. This contribution 
leads to a small shoulder as shown in Fig. ^ for £ = 0.6. 
For the larger elongation ( = 0.8, the shoulder gets more 
pronounced, and subsequently it leads to a broad peak 
at ri2 w 1.1 for C = 1-0. 

Unlike hard spheres for which two centers cannot ap- 
proach closer than r\2 = 1, two hard dumbbells can reach 
center-to-center distance ryi < 1 by adopting a "crossed" 
orientation. This explains why the g%${r\2) in Fig. ^ 
are positive even for ryi < 1- The probability of the 
"crossed" orientation increases with increasing density 
because it relieves the strain of the closely packed sys- 
tem. This effect is more pronounced at high elongations 
as shown in Fig. 0, and is a major factor contributing 
to the growth of the shoulder in g%§(? '12) for the region 
r\2 < 1 with increasing the elongation. 

Let us add one final comment. It is found from the 
extensions of the lower panel of Fig. |?] to larger 7*12 region 
that both the "T-shaped" and "CP-type" correlations are 
rather long-ranged. This is a manifestation of the strong 
peak in S^^q) at q i=s shown in Fig. || for ( = 1.0, as 
was discussed also in Refs. |l^ and g0[ The oscillatory 
feature of the above mentioned angular correlations are 
found to continue up to ri2 ~ 2ir/Aq, where Aq denotes 
the half width of that peak in S^iq)- This intermediate- 
range order is absent in the case of small elongations, say 
C = 0.4. 
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C. Bonding effects 

The dumbbell liquid for £ = 1.0 can be viewed as a 
system of hard spheres of diameter d — 1 whose density 
is 2p and where some additional covalent interaction has 
forced pairs to be formed. Let us consider the difference 
between the hard-sphere system and the bonded system 
in detail. Figure || compares the the radial-distribution 
functions and the static structure factors for the two sys- 
tems at the fixed packing fraction (p = 0.56. Notice 
that the site-site radial-distribution function g a b(f) for 
the symmetric dumbbell system becomes independent of 
the site indices a and b, and this is the adequate quan- 
tity to be compared with the radial-distribution function 
g(r) for the hard-sphere system. On the other hand, the 
total density static structure factor for the dumbbell 
system is the relevant one to be compared with the static 
structure factor S q for the hard-sphere system. This is 
because, when the packing fraction is fixed, the number 
density for the hard dumbbells with £ = 1.0 is half of 
that for the hard spheres, and the function properly 



accounts for this difference. The functions S q ~ ana s q are 
also the relevant inputs for the MCT equations for the 
hard-dumbbell and hard-sphere systems, respectively. 

It is seen from the upper panel of Fig. || that the agree- 
ment of the radial-distribution functions for the two sys- 
tems is very good except for the first-coordination-shell 
region. To demonstrate that this difference is primarily 
caused by the bonding, we shall consider the coordina- 
tion number K, i.e., the number of nearest neighboring 
spheres surrounding a central sphere, which can be cal- 
culated from the radial-distribution function. It is found 
for the hard-sphere system that K = 12.1 at ip = 0.56, 
which is a typical value for simple systems at high den- 
sity. So, K would tend to a value 12 also for the molecular 
system if the second sphere in one molecule was not at- 
tached to the first sphere. However, we found K = 11.4 
for £ = 1.0 at <p = 0.56. Thus, as should be expected, 
the second sphere in one molecule excludes one sphere in 
another from being a nearest neighbor to the first sphere, 
and this results in the reduction of the radial-distribution 
function in the first shell region as shown in the upper 
panel. 

The found feature for the radial-distribution functions 
also explains the reduction of the first peak height in the 
static structure factor for £ = 1.0 compared to that for 
C = 0.0, as exhibited in the lower panel of Fig. g, Let us 
consider what would happen to the static structure factor 
for the hard-sphere system when a short-ranged attrac- 
tive force is added. This problem has been discussed for 
a square- well system j$7|. As demonstrated there, the 
attraction causes bonding, in the sense that the most 
probable separation of two particles is smaller than ex- 
pected for a pure hard-sphere system. This leads to the 
shift of the first peak position in the static structure fac- 
tor to higher q, the decrease of the peak height, and the 
increase of the peak wings |5?J . Although the first fea- 
ture is not so prominent, the static structure factor for 
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FIG. 8: Upper panel: the site-site radial-distribution func- 
tion g a b(r) for the symmetric-hard-dumbbell system (see text) 
with ( = 1.0 at the packing fraction tp = 0.56 (solid line), 
and the radial-distribution function g(r) for the hard-sphere 
system (£ = 0.0) at the same packing fraction (dotted line). 
Lower panel: the corresponding static structure factor Sq 
for the total density fluctuations for the symmetric-hard- 
dumbbell system, and the static structure factor 5*, for the 
hard-sphere system. 



the £ = 1.0 dumbbell molecules reflects these features 
when compared to that for the hard-sphere system. 

We conclude that the structure of the cage for the hard- 
dumbbell system with £ = 1.0 is very close to the one for 
the hard-sphere system, and that the difference can be 
explained as being due to the bonding effect. 



IV. STRUCTURAL ARREST 



A. Critical nonergodicity parameters 

The upper panel of Fig. |^ exhibits the results for 
the nonergodicity parameters f q c at the critical point 
ip = (p c for the elongations £ — 0.4 and 1.0 calculated 
from Eq. (pi]). These are Debye- Waller factors of the 
system. They can be measured, in principle, as cross 
section for coherent neutron scattering. For large , 
the compressibility large. Therefore, sponta- 

neous arrest is easier for larger , and f q c exhibits a 
maximum near the first peak position of S q . With vary- 



ing q, fq oscillates in phase with S q (c/. Fig. ||). If the 
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FIG. 9: Critical nonergodicity parameter f q c and critical FIG. 10: Critical Lamb-Mossbauer factors for the tagged 



amplitude h q for the coherent total density fluctuations for 
the elongation £ = 1.0 (full lines), £ = 0.4 (dashed lines), and 
for the hard-sphere system (£ = 0.0, dotted lines). 



packing fraction increases, the arrested glass structure 
stiffens, i.e. the f{* increases. Expanding this increase 
for small distance parameters e = (ip—(p c ), one finds p, 0] 



f g N = ft + Dy/(<p <p c ) h N q + <p c ). (29) 

The critical amplitude is positive. It characterizes 
the susceptibility of the arrested structure with respect 
to changes of the control parameters. The formulas for 
the evaluation of and of the constant D > will be 
considered in the subsequent paper [p6| . Since < 1, 
fq ~ fq C i s bounded by 1 — f^ c . Therefore, the critical 
amplitude for the increase of is much smaller for 
q « 7 than for q off the structure-factor-peak position, 
as shown in the lower panel of Fig. ^. These features are 
analogous to those found in the hard-sphere system ||, |) . 

Figure flO| exhibits the tagged-molecule's critical noner- 
godicity parameters (x = N, Z) for the elongations 
C = 1-0 and 0.4 calculated from Eq. (|24|). These are 
Lamb-Mossbauer factors describing the arrested proba- 
bility distribution of the tagged molecule. As expected 
for a localized-distribution Fourier transform, the /" s - 
versus-g curves decrease with increasing q. The critical 
Lamb-Mossbauer factor for the total number density 
fluctuations approaches unity for q tending to zero due 
to the particle-number conservation law. On the other 
hand, there is no analogous conservation law for x = Z , 
and therefore one gets /£^ s < 1. 

A remarkable feature of Fig. [II] are the gentle oscilla- 
tions exhibited by f^f and f^ c s . In analogy to the dis- 



molecule's total density fluctuations f q J (upper panel) and 
for the charge-density fluctuations f q ^ (lower panel) for the 
elongation ( — 1.0 (full lines) and £ = 0.4 (dashed lines). The 
arrows indicate the wave numbers for the minimum positions 
in Wg and w q , respectively, calculated from Eq. (E6I). 



cussion for , it is expected that f^ c s oscillates in phase 
with Wg. As discussed in connection with Eq. (p6[), the 
function exhibits minima because of various ji(qC/2) 2 
contributions, and these minimum positions are marked 
as arrows in Fig. 10. One finds that the positions of the 



oscillations are well reproduced by the arrows, indicat- 
ing that they are due to the presence of various angular- 
momentum-index £ contributions to intramolecular inter- 
ference effects. A more definitive analysis concerning the 
origin of the oscillations should be based on the decom- 
position of fq C s in terms of the nonergodicity parameters 
of the tensorial density correlators <I>° u , (q, t) introduced 
in connection with Eq. fll2] ) p4| . Under the diagonal 
approximation &° sU ,(q,t) « 5u'&° s «(<?, t), one gets from 
the long-time limit of Eq. (021) : 



Jg,. 



(2/<)]T(2r 



l).M<zC/2) 2 / s c M,0). (30) 



Here, f^(q,£,0) = lim t ^oo $° w *) for V = fc and 
£ should be even (odd) for x = N (Z). As dis- 
cussed in Ref. Q the gentle oscillations exhibited by 
fq° can be explained as being due to interference ef- 
fects of the fg(q,£,0) with the intramolecular form fac- 
tors jg(q£/2) 2 . Unfortunately, we do not have informa- 
tion on the fg(q, £, 0). 

It might be interesting to consider how the results for 
the total density and charge density fluctuations can be 
translated to those for the atomic density fluctuations. 
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The latter can be obtained from the former via the in- 
verse relation of Eqs. ( fl6|) and (|l7|): 



R 



F AA = 



qN , rZ Z 



), 



N w N 

q,s q 



(31a) 
(31b) 



Notice that the self part F AA can be measured as cross 
section for incoherent neutron scattering. The results 
at the critical packing fraction are exhibited as solid 
lines in Fig. [ll] for the elongation £ = 1.0. The dashed 
and dashed-dotted lines are contributions from the total 
density (N) and charge density fluctuations (Z), respec- 
tively. A small peak centered at q « 4 develops in F AAc 
due to the charge-density fluctuations. The dotted line in 
the lower panel for the self part F AAc denotes the result 
based on the Gaussian approximation: 



F, 



AAc 



(32) 



Here, r A is the critical localization length for atom A, 
defined via lim g ^ (l — F AAc )/q 2 = (r A ) 2 . It is seen that 
the constituent atom's critical Lamb-Mossbauer factor 
F AAc is well described by a Gaussian, in particular, it 
does not exhibit oscillations. It is surprising that the 
sum of two non- Gaussian functions is almost Gaussian. 
The analogous results for £ = 0.4 are quite similar, except 



the peak of F AAc for q w 4 is suppressed. 
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FIG. 11: The full lines exhibit the critical nonergodicity 
parameters for the atomic density fluctuations F AAc (upper 
panel) and its self part F AAc (lower panel) for the elongation 
£ = 1.0. Fhe dashed and dashed-dotted lines denote con- 
tributions from the total-density (N) and charge-density (Z) 
components, respectively. The dotted line in the lower panel 
denotes the result based on the Gaussian approximation for 
F AAc , Eq. (§). 



B. Phase diagram 



The phase diagram in Fig. |l| can be understood as re- 
sult of the control-parameter dependence of the structure 
factors which were explained in Sec. HI. A prominent 
feature is the maximum of the <p c - versus-^ curve near 
C = 0.43. This is because of the two different mecha- 
nisms for the structural arrest, one dominating for small 
and the other one for large elongations. As discussed in 
connection with Fig. ^, the glass transition is driven by 
the first-peak region in the static structure factor Sg , 
irrespective of the elongation. For small elongations, the 
peak is primarily determined by the center-of-mass den- 
sity fluctuations, and its strength becomes weaker with 
increasing the elongation as explained in connection with 
Fig. [| and the upper panel of Fig. |5|. Therefore, a rela- 
tively higher packing fraction is required to get into the 
glassy phase if the elongation is increased, and this ex- 
plains the increase of the ifi c {() curve for small elonga- 
tions. On the other hand, for large elongations, the first 
peak in is mainly determined by the 1 = 2 angular 
correlation, and its magnitude gets larger with increasing 
the elongation as also explained in connection with Fig. 
and the lower panel of Fig. ^[ Thus, a relatively lower 
packing fraction is required for the glass formation as the 
elongation is increased, and this explains the decrease of 
the (fic(C) curve for large elongations. As a result of these 
two competing mechanisms for the glass formation, the 
transition line <^ C (C) exhibits a maximum. 



Another remarkable feature results from the structure- 
factor-peak reduction due to bonding which was ex- 
plained above in connection with Fig. ^| This reduction 
stabilizes the liquid phase. As a result, the critical pack- 
ing fraction for elongation £ = 1, <^ c (£ = 1) ~ 0.56, is 
larger than the one for the transition of the hard-sphere 
system, tp c (( = 0) » 0.53. Combined with the results 
discussed in the preceding paragraph, this implies that 
for all < C < 1 the critical packing fraction of the 
hard-dumbbell system is larger than that of the hard- 
sphere system; the liquid phase gets expanded due to the 
formation of molecules. The increased- free- volume phe- 
nomenon due to the bond formation is consistent with 
the result discussed for a square- well system p7| . 

There are two alternatives for the glassy states, phases 
II and III in Fig. |l|, with respect to the charge-density 
dynamics of the tagged molecule. Phase II deals with 
states for sufficiently small £. There is such small steric 
hindrance for a flip of the tagged molecule's axis between 
the two energetically equivalent positions e s and — e s 
that Eq. @ for x = Z yields f^ s = 0. The dynam- 
ics of the charge fluctuations is ergodic. In particular, 
the dipole correlator relaxes to zero: C\. s {t — > oo) = 
where C M (i) = (e s (t) • e s (0)) = <^L , S (*)' @. For suffi- 
ciently large C, on the other hand, the steric hindrance for 
dipole reorientations becomes so effective, that also the 
charge fluctuations behave nonergodically. In this case, 
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FIG. 12: Critical nonergodicity parameters f q ^ and the 
critical amplitudes h^ s for the tagged molecule's total den- 
sity fluctuations along the type-B-transition line in Fig. ^ 
parameterized by the elongation £ for the wave numbers 
q = 3.4(a), 7.0(6), 10.6(c), 14.2(d), and 17.4 (e). The criti- 
cal amplitudes k, for q = 14.2 and 17.4 have been omitted 
to avoid the overlapping of the curves. 



FIG. 13: Results as in Fig. hi but for the tagged molecule's 
charge density fluctuations, j£ t s an( A hq,s- The dashed lines 
are added here, denoting the critical nonergodicity parameter 
(upper panel) and the critical amplitude (lower panel) for 
the zero- wave- number limit. The arrow marks the transition 
point from phase II to phase III at f c = 0.345 taken from 
Fig. g 



Eq. ( p4| ) for x = Z yields a positive long-time limit, < 
fq a = 4>q s (t — > oo). In particular, dipole disturbances 
do not relax to zero: Ci :S (t — > oo) = fi iS = fg— 0tS > 0. 
This phase III is a glass with all structural disturbances 
exhibiting nonergodic motion. In particular, the noner- 
godicity parameter /f s for, say, £ > 0.6 is as large as the 
maximum of f^ c , Fig. ||. The two phases II and III are 
separated by a curve <^a(C)j where <pa(C) > Vc(C)- This 
curve is shown as the dashed line in Fig. [j]. Since the 
steric hindrance for the molecule's flip motion increases 
with increasing £, one might expect that the curve <pa(C) 
would monotonically decrease with increasing £. How- 
ever, this is not the case. What monotonically decreases 
with increasing £ is the difference <Pa{C) ~ Vc(C)j an d 
the variation of <p c (C) dominates that of <pa(C) f° r small 
£ c — £. The curve ^a(C) terminates at the critical elon- 
gation ( c : <pa((c) — <Pc(Cc)- For our model, one finds 
Cc = 0.345, and its position is marked by an arrow in 
Fig. [lj The asymptotic laws for the transition from phase 
II to phase III have earlier been described as the type- 
A transition, as can be inferred from Ref. [38| and the 
papers quoted there. The square-root singularity of the 
Debye- Waller factor f*, Eq. @, implies via Eq. (g) 
that the two phase transition lines do not merge transver- 
sally: J=</?a(Cc) = 7%<Pc{(c)- All together, this explains 
the minimum of the <p.a(£)— versus— £ curve near £ = 0.23. 
There are some characteristic features of the type- A 



transition which are relevant in the analysis of the type- 
£?-transition dynamics: these are connected with the Co- 
variation of the critical Lamb-Mossbauer factors f^° s and 
the critical amplitudes h x q a [^6|. Figure [l2| exhibits 
and h^ s for the total density fluctuations along the type- 
B transition line <fi c {C) parameterized by £. Figure |i"3| 
shows the corresponding results for the charge-density 
correlator, / ? , and s . They deal with the transition 
from phase II for £ < £ c to phase III for £ > £ c . We note 
in passing that the curves shown in these figures exhibit 
non-monotonic £ dependence, such as wiggles, or even 
minima and maxima. These anomalies are analogues 
of the gentle oscillations, discussed above in connection 
with Fig. [l^, and they can in principle be explained as 
was done in Ref. |24| for dumbbell molecules immersed in 
a hard-sphere system. For strong steric hindrance, say 
C > 0.8, fgg is rather close to f^°, and this holds also 
for critical amplitudes, h qs and hg S . For £ approach- 
ing C c , the Lamb-Mossbauer factor falls below /^f, 
and the critical amplitude hg S grows above h^ s . These 
are characteristic features of the type-^4 transition, whose 
transition point can be characterized by the vanishing of 
the critical nonergodicity parameter j^ c s and by the di- 
vergence of the critical amplitude h^ s Q. The former 
feature is demonstrated, for example, by the strong de- 
crease of /f s — fq^Q s for C = 0.4 shown in the lower 
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panel of Fig. [H] compared to the one for ( — 1.0 in the 
same panel. Since the critical amplitude h^ s gauges the 
dynamics in the /3-relaxation regime, the dynamics of the 
charge-density correlators as well as the dipole correlator 
near Cc is strongly influenced by precursor phenomena of 
the type- A transition from phase II to phase III. Their 
dynamics in the a-relaxation regime is also perturbed 
since the leading correction to the a-scaling law is pro- 
portional to the critical amplitude These features 
will be discussed in the following paper pq| . Thus, the 
dynamics for elongations close to Cc is qualitatively dif- 
ferent from the one for large elongations, and this is why 
we have chosen as the representative elongation C = 0.4 
for the demonstration of the state with weak steric hin- 
drance for reorientations. Generically, there is no type- A 
transition line for arbitrary diatomic molecules. In our 
problem, this singularity is due to an additional sym- 
metry that produces vanishing coupling constants. The 
top-down symmetry of the molecule renders the MCT 
equations decouple completely into a set for total density 
fluctuations and another one for charge-density fluctua- 
tions. For nearly top-down symmetrical molecules, the 
type-j4 transition is smeared to a rapid crossover from 
the very small nonergodicity parameters f j?. for C <C Cc 
to such of order unity for ( ^> £ c ]2l], |8[ [59[. By conti- 
nuity, the critical amplitude s remains large for C near 
Cc also for asymmetrical molecules. Therefore, the re- 
sults for ( ~ 0.4 are also representative for such cases, 
where the type-A transition singularity is avoided due to 
breaking of the top-down symmetry of the constituent 
molecules, provided the breaking is sufficiently weak. 

The described type-A transition has been studied also 
for a single dumbbell immersed in a system of hard 
spheres. The diagrams corresponding to the up per pan- 
els in Figs. |l2| and [l3] are qualitatively similar jgjfl , but 
there are two remarkable differences. First, the form fac- 
tors f^g are somewhat larger and the variation with C 
for C > C c is more pronounced for the hard-dumbbell 
system than for the corresponding quantities for the sim- 
ple system. Similarly, for C > 0.5 the /^r are larger in 
Fig [[3] than for the single-dumbbell system. Second, the 
(C — Cc) interval for the decay of /^r from large weakly 
C-dependent values to zero at the transition is narrower 
for the motion in the dumbbell liquid than for the motion 
in the hard-sphere system. These differences reflect the 
fact that steric hindrance for translation as well as for 
reorientation is more efficient if the cage-forming neigh- 
bor molecules are sufficiently elongated rather than being 
spherical. This conclusion explains also that the found 
critical value Cc = 0.345 is smaller than the correspond- 
ing value 0.380 obtained for a single dumbbell in a hard- 
sphere system p4|. 



V. CONCLUSIONS 

A mode-coupling theory (MCT) for the evolution of 
jlassy dynamics is derived and used to discuss the ide- 



alized liquid-glass transition in a hard-dumbbell system 
(HDS). The theory predicts a singular change of the dy- 
namics caused by a regular change of the canonically de- 
fined equilibrium structure factors with variations of con- 
trol parameters like the packing fraction ip. The struc- 
ture factors define the mode-coupling constants in the 
equations of motion for the correlation functions, and 
they have been evaluated within the RISM and Percus- 
Yevick theories. The good agreement of the results from 
the two approximate approaches support the opinion 
that the used input information of the MCT is semi- 
quantitatively correct. The results have been used to 
demonstrate that "T-shaped" configurations are the pre- 
ferred arrangements of the cage-forming neighbors of a 
molecule. The found arrangements are similar to those 
discussed earlier for more dilute systems [ j33| ; but the or- 
dering in our high-density regime is more pronounced. In 
addition, there is intermediate-ranged order leading to a 
central peak for quadrupole-density fluctuations, but this 
is irrelevant for the explanati on o f the glassy dynamics 
within the present theory (Sec. III). There is no informa- 
tion available on the correctness of the cited structure- 
factor theories within the large-tp regime studied in this 
paper. This implies obvious reservation concerning quan- 
titative details of the results presented. 

The MCT for molecular systems proposed in this pa- 
per is based on describing the dynamics by n-by-n- 
matrix correlators formed with the n interaction-site den- 
sities of the molecule's constituent atoms. Such basis 
is inferior to the one using a description by infinite- 
matrix correlators formed with tensor-density fluctua- 
tions @, |l|, [l6[ H, ||, |(], || H , provided the equa- 
tions of motion of the latter theory could be solved for pa- 
rameters and time regimes of interest. For example, our 
theory does not directly lead to results for the angular- 
momentum-l = 2 reorientational correlator, which is rel- 
evant for the description of depolarized light-scattering 
data. However, it was shown already in some other con- 
text [^4|, how the £ = 2 reorientational correlator can 
be obtained as an addendum of the site-representation 
theory. A more subtle extension of the theory would be 
necessary, if there would be a second-order phase tran- 
sition. A treatment of the interference of the slow glass 
dynamics with the critical dynamics of the phase tran- 
sition would require the inclusion of the critical fluctua- 
tions in the relaxation kernel in the spirit of the original 
derivation of mode-coupling theories by Kawasaki [|40"|| . 
It is unclear at present whether such an extension can be 
formulated. 

Compared to a hard-sphere system (HSS), the fu- 
sion of two hard spheres to a dumbbell of elongation 
C, < C < I, increases the free volume if the packing 
fraction is kept fixed. Therefore, the liquid gets stabi- 
lized and the line for the liquid-to-glass transition ip c (C) 
is above the transition value tp^ ss of the HSS. Like for 
the HSS, the transition is driven by the density fluctu- 
ations with wave vectors near the position of the first 
sharp diffraction peak. For small C; the peak is formed 
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by the center-to-center correlations which decrease with 
increasing leading to an increase of the t/? c -versus-C 
transition curve. For large £, the peak is formed by the 
quadrupole correlations; and these increase with £, lead- 
ing to a decrease of the transition line. This explains the 
pronounced maximum of the transition curve in Fig. [l]. 
The model studied exhibits a symmetry with respect to 
the top-down flip of the molecule's axis. As explained 
in the earlier MCT literature, this implies a line of spin- 
glass-type transitions shown as dashed curve in Fig. [y. 

A comment concerning the accuracy of the reported 
calculations might be adequate. After the specified 
discretization of the wave numbers, Eq. (|2l| ) for the 

100 numbers f q is solved by the iteration = 
jrN [f NU)]/ {1 + F»[f»U)]}, j = 0,1,..., starting 

from fq = 1. The sequence decreases monotoni- 
cally towards the nonergodicity parameter fq^ — > f¥. 
The linearized iteration for Sfq — fq''' — reads 
5fq +1 ^ = ~}2 p Aq P 8fp \ where the Frobenius matrix A 
is given by A qp = (1 - f» ) 2 dF™[f N ]/df£ . The matrix 
has a maximum eigenvalue E < 1. Off the critical points, 
one gets E < 1, and the convergence of the iteration is 
exponentially fast. The critical point is characterized by 
E c = 1, and here the convergence is only algebraically. 
The proofs of the cited mathematical properties can be 
found in Ref. [oj. Near the critical point, one derives from 
Eq. (|29|): E c — Eoc y/ip — tp c . Analogous statements hold 
for the calculations of fg S from Eq. ( pif ) . In our numer- 
ical work, the value for E is controlled and ip — tp c is 
determined routinely so, that E c — E w 10~ 4 . Hence, 
the critical points are calculated with an accuracy of the 
order 10 -8 . Thus, the accuracy of the lines in Fig. [I] is 
determined by the number n* of values for £ used to cal- 
culate f c (C) an d <Pa(0- We used a grid with n* = 70; 
it was chosen non-uniformly over the interval < Q < 1 
with the highest density of points for £ near ( c . 

Letz et at have discussed a liquid-glass phase diagram 
for a system of hard ellipsoids |19| . Considering their as- 
pect ratio of the prolate ellipsoids as analogue of (Q + 1) 
for the dumbbells, their phase diagram looks similar to 
Fig. |l|. They also show the analogue of the glass-to-glass 
transition curve, albeit without a minimum and with a 
transversal termination at the liquid-glass-transition line. 
It is argued in Ref. |l^ that the strong decrease of the p c ~ 
versus-^ curve for aspect ratios near 2 is an implication of 
the central peak of the angular-momentum-^ = 2 corre- 
lations, reflecting a nematic-transition precursor. Hence, 
the explanation of the phase diagram given in Ref. [l^ 
for the hard-ellipsoid system is quite different from the 
explanation of Fig. [j] for the dumbbell system. 

The critical form factors for the glass fg quantify 
the arrested amorphous density fluctuations at the tran- 
sition. The wave-vector dependence, Fig. 0, is quite sim- 
ilar to that for a HSS. This reflects the fact that the 
cage around an interaction site is quite similar to the one 
found for the HSS. The critical nonergodicity parameters 



fqg for the arrest of the dipole reorientation, Fig. |l^, are 
larger than the same quantities calculated for a single 
molecule in the HSS Q. In particular, the decrease of 
fq=o s f° r C decreasing to the critical value £ c is so abrupt, 
that the transition looks similar to a discontinuous one. 
This shows that steric hindrance for reorientation is more 
effective in a molecular system than in a system of spher- 
ical particles. 

In the following paper (^6), it will be shown that the 
results for the arrested structure provide the key for an 
explanation of the structural relaxation. 
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APPENDIX A: MODE-COUPLING THEORY FOR 
MOLECULAR SYSTEMS 

The MCT focuses on the dynamics of density fluctua- 
tions. Within the site representation, the basic variables 
are the density fluctuations for the n interaction sites of 
the N molecules: pi — J2iLi exp(i<f • r* a ), a = 1, 2, • • • , n. 
Here r° denotes the position vector of the site a in 
molecule i. The most important correlation functions 
for a statistical description of the dynamics are 

Ff{t) - (p|(t)V|(0))/Ar, a, b = 1, ■ ■ ■ ,n. (Al) 

These n 2 functions shall be considered as the elements of 
an n x n matrix F q (t) . This matrix is real and symmetric. 
The short-time expansion of this matrix is given by 

¥q{t) = S q -\ q 2 3 q t 2 + 0(t 4 ). (A2) 

Here the structure-factor matrix is given by S® — 

(Pm*p}(0))/N, and Jf - ((q ■ j q ?y(q- j*))/Nf is 
defined in terms of the currents referring to the inter- 
action sites jg = YliLi exp(iq ■ r") with vf denoting 
the velocity of the site a in molecule i. The Zwanzig-Mori 
formalism [ p8| leads to an exact equation of motion for 
F q (t): 

dfF q (t) +n 2 q F q(t)+n 2 q f dt , m q (t-t')d t ,Fq(t')=0, 

Jo 

(A3) 

where 

n 2 l =q 2 3qS-\ (A4) 

The right-hand side of this equation is a product of two 
positive definite matrices. Hence it is equivalent to the 
square of a positive definite matrix. Therefore, one can 
write it as the square, ft 2 , of some matrix Q q . Splitting 
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off this matrix in front of the convolution integral is done 
for later convenience. 

The difficult problem is the derivation of an approx- 
imation for the matrix m q (t) of fluctuating- force corre- 
lations such that the cage effect is treated reasonably. 
This has been done originally in Ref. |23l by extending 
the procedure used for atomic systems 142], But the re- 
ported formulas |^3| are not acceptable. Firstly, they do 
not properly reduce to the ones for simple systems in the 
united atom limit. Secondly, the momentum conserva- 
tion law for coherent dynamics is not satisfied. For these 
reasons, an alternative derivation has been developed in 
Rcf. ^4] starting from the projection-operator theory of 
Mori and Fujisaka [^3|, albeit for molecules immersed in 
a simple system. It is possible to generalize this deriva- 
tion for the coherent density correlators F q (t) for molec- 
ular systems, but the procedure becomes more involved; 
it shall be described in a separate paper Q|. Here, a 
more simplified derivation shall be presented. 

The simplified procedure starts by assuming that a 
molecular system is a mixture of constituent atoms; in- 
tramolecular constraints between constituent atoms are 
accounted for by the pair correlations only. In this way, 
a system of iV molecules is treated as a mixture of n 
species, each consisting of N particles. Using the equa- 
tions in MCT for mixtures ]45[ ], one gets the following 
expression for the relaxation kernel: 



mf{t)=Tf [¥{£)], 



(A5) 



where the mode-coupling functional J- q is given by the 
equilibrium quantities: 



x/j^'ir'. ( Ae ) 

<w(<7l k,p) = {q- [<*"" fccf + S aX pc a /}} 

x {q- [6 b »' kef + 5 bX ' pc b /}} I q\ (A7) 

with p = q — k. Here, the direct correlation function 
is defined via the Ornstein-Zernike equation for a mix- 
ture ||, pcf = S ab - [S q ~ 1 ] ab . Now, let us turn on the 
intramolecular constraints between constituent atoms. 
This amounts to replacing the direct correlation func- 
tion c° h for a mixture of spherical particles with the one 
for molecular systems defined via the site-site Ornstein- 



Zernike equation 127 



ab 



-lyab 



is, 



-llafc 



,,a& 



Here 
de- 



enter the intramolecular correlation functions 
scribing the constraints. The so obtained equations for 
a mixture contain a frequency matrix Q, 2 which reflects 
the 3n independent degrees of freedom of the molecule. 
In particular, fl 2 gets an n-fold degenerate eigenvalue 
zero for q = due to the n particle-number-conservation 
laws for the n species. The used classical theory cannot 
account for the fact that vibrational degrees of freedom 
are frozen out at sufficiently low temperature because of 



quantum effects. To repair this shortcoming, we make 
the assumption that in the regime of interest the rigid- 
ity of the molecule can be accounted for by replacing the 
classical flexible-molecule value for £l 2 in front of the con- 



volution integral in Eq. (A3) by the formula in Eq. (A4). 



The matrix fl 2 for a rigid molecule exhibits only one 
eigenvalue zero for q = 0, since there is only one inde- 
pendent conservation law for the number of molecules. 

The MCT equations for the tagged-molecule density 
correlator F qs (t) defined in Eq. (||) can be obtained sim- 
ilarly, and only the resulting equations shall be quoted. 
The exact Zwanzig-Mori equation reads 

d?F q , s (t) + nl s F q , s (t) 

+S1 2 I dt'm q , s {t-t')d t >F q , s (t') = Q, (A8) 
Jo 



whe re th e characteristic frequency matrix is given as in 
Eq. ( A4) by tt 2 s = q 2 3 q w" 1 . The expression for the re- 
laxation kernel can be formulated as mode-coupling func- 
tional !F a s : 



mt(t)=^\[F s (t),F(t)). 



(A9) 



The explicit expression for the functional F q%s reads, 
with p = q — k. 



F q h s\fsJ\ 



(2tt 



,3 



ac 

IT 

X.jj, 



dk 



q-p 



xfgjp. (A10) 



Let us note some mathematical results valid for the 
MCT formulated above. First of all, there is a solution 
F q (t) of the nonlinear equations of motion for all times 
t. This solution is uniquely fixed by the initial condi- 
tions F q (t = 0) = S q and d t F q (t = 0) = 0. For every 
finite time interval, the solution depends smoothly on the 
numbers tt 2 ab and V£ b , , (q; k, p) . The solutions are cor- 
relation functions in the sense that they can be Laplace 
transformed to functions having a spectral representation 
with a spectrum F^'(w) which is a positive definite ma- 
trix. The matrix of long-time limits F q = lim^oo F q (t) 
obeys the set of implicit equations defined by the mode- 
coupling functional !F q : 



[S, - F q ] 



(All) 



Let us remember that there is a semi-ordering in the 
space of real symmetric n-by-n matrices, A > B, de- 
fined by A — B to be positive definite. With this no- 
tation, the maximum theorem holds: if F q is a solu- 

= T q [F], then 



tion of Eq. (All), i.e. if F q [S q - F q \ 



F q < F q . If an iteration sequence F q ' , j = 0, 1, ■ • • , is de- 



fined by F 



O+i) r 



P (0) 



S q , then F q 



(j+i) 



g 

C Fjj^ and lim 



J- q \F^\ starting from 



?0) 



If the Jacobian of Eq. (All) does not have a vanish- 
ing eigenvalue, the long-time limit F q depends smoothly 
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on the coupling coefficients V£y ,(q; k,p). A singular- 
ity can occur only if an eigenvalue vanishes. The sub- 
tlest property is that such vanishing eigenvalue is non- 
degenerate. Hence, using the terminology of Arnold jh|, 
all possible singularities are bifurcations of the cuspoid 
type Ai, £ = 2, 3, • • • . The generic singularity for changes 
of a single control parameter is, as for the MCT of simple 
systems, a fold bifurcation A 2 . It is then obvious, that all 



universal results for simple systems are valid also for the 
MCT for molecular systems formulated above. For exam- 
ple, Eq. ( p9| ) holds with and replaced by positive 
definite matrices. The proofs of the cited results of this 
paragraph shall not be described here for brevity, since 
essentially the same issues for matrices of density corre- 
lators have been independently discussed and proved by 
Franosch and Voigtmann f47|| . 
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